Light cluster production in E/A = 61 MeV 36 Ar + 112124 Sn reactions 
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Experimental kinetic energy distributions and small-angle 
two-particle correlation functions involving deuterons and tri- 
tons are compared for 36 Ar+ 112,124 Sn collisions at E/A = 
61 MeV (i.e. for systems similar in size, but with different 
isospin content). A larger triton yield is observed from the 
more neutron-rich system, as predicted by IBUU simulations, 
while the emission times of the light clusters are found to be 
the same for the two Sn-target systems. For both systems, 
the time sequence r<j < r p < it, is deduced for charged parti- 
cles emitted from the intermediate velocity source. 
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I. INTRODUCTION 

Among the challenges in nuclear physics today, is to 
understand the nuclear interaction for systems with an 
exotic composition of neutrons and protons (which may 
occur in the r-process of nucleosynthesis and in neutron 
stars). One of the topics in this challenge is to understand 
the isospin dependence of the nuclear equation of state 
(EOS), and in particular its density dependence [1-6]. In 
this difficult task many observables are needed to put as 
many constraints as possible on the nuclear interaction. 
Among several suggested observables, [7-13], the two- 
nucleon correlation function has been considered in Rcfs. 
[14,15]. 

Experimentally, isospin effects in two-particle correla- 
tion functions were investigated by our group in Ref. [16]. 
The experimental data were from E/A = 61 MeV 36 Ar- 
induced semi-peripheral collisions on isotope-separated 
H2,i24g n targets. Stronger np correlations were found 
for the Ar + 124 Sn system, pointing to a shorter neu- 
tron emission time from the more neutron-rich system. 
Smaller isospin effects were seen also in pp, nd and nt cor- 
relation functions [16]. In addition, a study of the emis- 
sion time sequence of neutrons, protons, and deuterons 
from the same reactions was carried out by our group in 
Ref. [17], utilizing the method of the velocity-gated cor- 
relation functions. That investigation revealed a sensitiv- 



ity of the particle emission time sequence to the isospin 
content of the emitting source, particularly for intermedi- 
ate velocity source emission and emission from the target 
residues. 

In this paper we present complementary information 
to that of Refs. [16,17], from the same experimental data 
set, focusing on deuterons and tritons. Indeed, also light 
cluster production is expected to be a sensitive probe 
for the isospin dependence of the EOS. In the theoreti- 
cal investigation of Ref. [18] (where IBUU was coupled 
with a coalescence model for cluster production), it was 
found that the nuclear symmetry energy significantly af- 
fects the production of light clusters in heavy ion col- 
lisions. More deuterons and tritons are produced with 
a stiff nuclear symmetry energy than with a soft. This 
is because the stiff symmetry energy induces a stronger 
pressure in the reaction system, and thus causes an ear- 
lier emission of neutrons and protons leading to stronger 
correlations among nucleons, and, consequently, larger 
yields of deuterons and tritons. 



II. EXPERIMENTAL DETAILS 

The experiment was performed at KVI (Groningen). 
The interferometer consisted of 16 CsI(Tl) detectors for 
light charged particles, mounted in the angular range 30° 
< < 114°, and 32 liquid scintillator neutron detectors, 
mounted behind the "holes" of the Csl array, in match- 
ing positions to provide the np interferometer [19]. En- 
ergy thresholds for protons, deuterons, and tritons in the 
CsI(Tl) detectors were 8, 11, and 14 MeV respectively. 
Details about the experimental setup and the particle 
energy determination are given in Refs. [19-21]. 

From our previous analyses [16,17], including the 
source analysis of Ref. [21], we have confirmed that 
the emission of light particles from 61 A MeV (semi- 
peripheral) heavy ion reactions originates from (at least) 
three sources: projectile and target residue evaporative 
sources (statistical evaporation) and intermediate veloc- 
ity source. The intermediate source represents dynami- 
cal emission, which is described either by early nucleon- 
nucleon collisions or by other pre-equilibrium processes 
[22-30] . Within the multi-source reaction mechanism de- 
scribed above, composite particles, like deuterons and tri- 
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tons, are believed to be predominantly emitted from the 
dynamical emission source [31,32], where they are formed 
by a coalescence mechanism [21,33]. 



III. RESULTS 

Figure 1 compares the kinetic energy distributions 
of deuterons and tritons measured in 36 Ar+ 112 Sn and 
36 Ar+ 124 Sn collisions. One can notice that the yields 
of the deuteron spectra are the same for the two Sn- 
targets, while the triton yield from the more neutron-rich 
36 Ar+ 124 Sn system is enhanced (by a factor 1.1 to 1.5) 
over the whole angular and energy range considered. 
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FIG. 1. Comparison of the d and t kinetic energy yields 
measured in 36 Ar+ 124 Sn (dashed lines) and 36 Ar+ 112 Sn (full 
lines) at the angles indicated in the figure. The yields are 
arbitrarily shifted in the y-axis as indicated in the figure. 

To further investigate the deuteron and triton produc- 
tion, we present the two-particle correlation functions in- 
volving deuterons and tritons. The correlation function, 
C(q, Ptot) = kN c (q,P tot ) / N nc (q,P tot ), is generated by 
dividing the coincidence yield (iV c ) by the yield of non- 
correlated events (N nc ) constructed from the product of 
the singles distributions [34]. The relative momentum is, 
q = fi(pi/mi — ^2/^2), where a is the reduced mass, 
and the total momentum is, P tot = Pi + P2- The corre- 
lation function is normalized to unity at large values of 
q, (200 < q < 300 MeV/c for all particle pairs) where no 
correlations are expected. 

The pp, np, and pd correlation functions have already 
been presented in Refs. [16,17]. In those analyses, par- 
ticles emitted by the intermediate velocity source were 
enhanced by selecting high- and intermediate- P to t pairs, 
while particles emitted by the target residues were en- 
hanced by selecting low-P tot pairs in the target frame. 
Further selection of different sources was obtained by ap- 



plying angular gates. The pp and np correlation functions 
were found to be very sensitive to the applied gates, indi- 
cating that both neutrons and protons are emitted from 
several sources. In contrast, deuterons and tritons mainly 
originate from common sources with a small spread in 
source velocity (v source ~ 0.18 c) [21]. In agreement with 
those findings, we find, in the present work, that the dd, 
tt, and dt correlation functions are quite insensitive to 
applying total-momentum and angular gates. We have 
therefore chosen to present angle and total-momcntum- 
integrated dd, tt and dt correlation functions in this pa- 
per. 

Figure 2(a,b) compares the dd (a) and tt (b) correla- 
tion functions measured from the two Sn-target systems. 
One can see that the Coulomb hole dominates both corre- 
lations at low values of the relative momentum. For both 
dd and tt, the width of the Coulomb hole is the same for 
the two Sn-targets. This feature indicates that the av- 
erage emission time of deuterons and tritons is the same 
for the two systems. The width of the Coulomb hole is 
larger for dd than for tt, which is consistent with a larger 
average emission time for tritons than for deuterons. 
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FIG. 2. From 36 Ar+ 124 Sn (open circles) and a6 Ar+ 112 Sn 
(filled circles), dd (a), tt (b), and dt (d) angle- 
and total-momentum-integrated correlation functions. A 
high-total-momentum gate is applied to the pt (c) correlation 
function. 

The dt correlation function in Fig. 2d is very similar for 
the two Sn-targets. It is governed by strong resonances at 
q « 10.8 MeV/c and q « 83.5 MeV/c, due to the excited 
states, E* = 16.75 MeV (T cm = 0.076 MeV) and E* = 
19.8 MeV (r cm = 2.5 MeV) of 5 Hc [35]. 

We now turn to the correlation functions involving 
also protons. The pd correlation function, that was pre- 
sented in Refs. [16,17], is characterized by a pronounced 
anti-correlation at small relative momenta, due to fi- 
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nal state Coulomb interaction. The anti-correlation is 
more pronounced for particle pairs selected by high- and 
intcrmediate-total-momentum gates, than for pairs se- 
lected by low-total-momentum gate. 

The pt correlation function, shown in Fig. 2c, contains 
contributions of excited states of 4 He at q « 23.6 MeV/c 
(E* = 20.21 MeV, Y cm = 0.5 MeV), q « 41.0 MeV/c (E* 
= 21.01 MeV, T cm = 0.84 MeV), and q w 53.4 MeV/c 
(£* = 21.84 MeV, T cm = 2.01 MeV) [36], while higher 
excited states of 4 He are much broader and not visible in 
Fig. 2c. The positions of the 1 st and 3 rd resonances are 
marked by arrows. 

In contrast with dd, tt, and dt, the pd and pt correla- 
tion functions are sensitive to applying angle- and total- 
momentum gates. This is because protons are emitted 
from several sources [21]. In order to consistently inves- 
tigate emission from the intermediate velocity source, we 
have therefore applied a high-total-momentum gate to 
the pt correlation function (shown in Fig. 2c). This se- 
lection generates a correlation function characterized by 
stronger resonance peaks as compared to those obtained 
by lower total momentum gates. Inspection of Fig. 2c re- 
veals also a small difference in the correlation functions 
for the two Sn isotopes. The more neutron- rich system 
leads to a slightly stronger resonance peak for the higher 
excited states of 4 He, while the opposite behavior is seen 
for the lowest excited state. 

Further information about the emission times and 
chronologies can be obtained by means of particle- 
velocity-gated correlation functions of non-identical par- 
ticles [37-39]. This adds an important and novel piece 
of information to the picture of the reaction mechanism 
emerging from previous studies [16,17]. For non-identical 
particles, a and b, we construct the correlation functions 
C a (q), gated on pairs with v a > Vb, and Cb(q), gated 
on pairs with Vb > v a . The particle velocities are calcu- 
lated in the frame of the emitting source, and the same 
normalization constant, calculated from the ungated cor- 
relation function, is utilized for both C a and Cb [38]. If 
two particles a and b are emitted independently from a 
source, with a later (earlier) than b, then the ratio C a /Cb 
will show a peak (dip) in the region of q where there is 
a correlation, and a dip (peak) in the region of q where 
there is an anti-correlation. 

To correctly interpret the experimental results, it is im- 
portant to realize that the resonance peaks may a priori 
have two different origins [40,41]: 

1 . From processes where an unstable fragment formed 
in the reaction decays into two particles which are 
measured in coincidence (e.g. 5 Hc — > d + t). The 
energy and momentum conservation in the decay 
determines the location of the resonance peak. 

2. From final state interactions between particles 
emitted independently from a source (e.g. in- 
termediate velocity source, target and projectile 
residues). 



If the two particles are emitted independently, the 
velocity-gated correlation functions contain information 
on the time sequence, but this is not the case if the two 
particles originate from the two-body decay of an un- 
stable fragment emitted in the reaction. In the latter 
case, the two particle velocities are determined by mo- 
mentum conservation, and in the rest system of the de- 
caying fragment, the lightest particle will always get the 
highest velocity. In this case, the velocity-gated correla- 
tion function obtained with the condition that the light- 
est particle has the largest velocity, should exhibit the 
strongest correlation or anti-correlation (see Ref. [42] for 
further details). When we in some cases observe the op- 
posite behavior (namely that the gate, where the heaviest 
particle has the largest velocity, leads to a stronger corre- 
lation or anti-correlation) we can reliably conclude that 
this behavior is dominated by a mechanism other than 
two-body decay. We attribute this effect to the interac- 
tion of independently emitted particles, and in this case, 
we use the velocity-gated correlation function to obtain 
information on the time sequence of the independently 
emitted particles. 
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FIG. 3. From E/A = 61 MeV 36 Ar + 112 Sn (loft column) 
and 124 Sn (right column) collisions, particle-velocity-gated 
(filled and open circles) pt (a,c) and dt (e,g) correlation func- 
tions, and their ratios (b,d,f,h). 

The study of the particle- velocity-gated pd correlation 
function can be found in Ref. [17]. The extracted emis- 
sion time sequence is that deuterons are, on average, 
emitted earlier than protons for both Sn-target systems 
(consistent with the findings of previous experimental 
investigations [38,39]). Note that this result refers to 
the average emission time, and it does not exclude that 
prompt protons are emitted before any deuteron is emit- 
ted. 

The particle-velocity-gated pt correlation function is 
shown in Fig. 3a-d. In the momentum region below q w 
50 MeV/c, the enhancement of C p (v p > v t ) over C\ (i>t > 
v p ), expected from momentum conservation of the two- 
body decay of 4 He, is not observed. This indicates that 
the behavior of the pt correlation function in this region is 
dominated by independently emitted pairs. The emission 
chronology that can be deduced from this region is that 
protons are emitted earlier than tritons. Inspection of 
the Cp/Ct ratios (Fig. 3b, d) reveals that the deduced pt 
emission chronology is qualitatively similar for the two 
Sn-target systems. 

The particle-velocity-gated dt correlation function in 
Fig. 3e-h, shows an enhancement of the Cd correlation 
function that is expected from the two-body decays of 
5 He. The behavior of the Cd/C t ratio (Fig. 3f,h) is there- 
fore of delicate interpretation, since it may be dominated 
by the two-body resonance decays of 5 He, and it does not 
led itself to deduce the emission chronology in this case. 



IV. SUMMARY 

Kinetic energy spectra and two-particle correlation 
functions of deuterons and tritons have been investigated 
in E/A = 61 MeV 36 Ar+ 112 < 124 Sn reactions. The yield 
of tritons is larger from the more neutron-rich system, 
while no differences are found in the deuteron yield. The 
dd, tt and dt correlation functions are rather insensitive 
to applying angle- and total-momentum gates, support- 
ing the interpretation [21] that the origin of emission of 
deuterons and tritons is a coalescence mechanism tak- 
ing place in the intermediate velocity source created in 
the early dynamical phase. For pd and pt pairs, instead, 
the emission from the intermediate velocity source can 
be enhanced by selecting high-total-momcntum pairs. 

From Ref. [17] and the analysis of the velocity- gated 
pt correlation function presented here, we find that, for 
the intermediate velocity source, deuterons are on aver- 
age emitted before protons, which in turn are emitted 
before tritons, Td < r p < r t . This is also consistent with 
the behavior of the dd and tt correlation functions. 

The isospin effects in deuteron and triton emission, 
appear to be weak for the studied systems. Apart for 



a larger triton yield from the more neutron-rich system, 
the emission chronology appears to be the same, as well 
as the average deuteron and triton emission times. One 
possible interpretation is that the larger yield of tritons 
in the 36 Ar+ 124 Sn system is formed in the same amount 
of time as in the 36 Ar+ 112 Sn system, because of a larger 
abundance of neutrons emitted early in the 36 Ar+ 124 Sn 
reaction. 
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